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Abstract 

We study massless g^^^-theory at high temperature and with zero chemical poten- 
tial. Using modern effective field theory methods, we derive the screening mass squared 
to order g"^ . It is demonstrated that the result is rcnormalization group invariant to 
this order in the coupling constant. Rcnormalization group methods are used to sum 
up leading logarithms of g. 

PACS number(s): ll.lO.Wx, 12.20.DS, 12.38.Bx 



1 Introduction 

There has been tremendous progress in perturbative calculations in thermal field 
theory, since the work on symmetry behaviour at finite temperature by Dolan and 
Jackiw more than two decades ago 0]. 

One major breakthrough was the discovery by Braaten and Pisarski that ordi- 
nary perturbation theory breaks down high temperature for soft external momenta 
(fc ~ gT, where g is the coupling constant) [Q. The usual connection between 
the number of loops and powers of the coupling constant is lost. The leading order 
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results for physical quantities receive contribution from all orders in perturbation the- 
ory. The solution to the problem is to use an effective expansion which resums this 
infinite subset of diagrams and where loop corrections are truly perturbative (down 
by powers of g). This resummation program and extensions thereof, have been ap- 
plied to calculate a variety of static and dynamical quantities in hot plasmas (see 
e.g. Ref. [3,4]). The literature is now vast, but we shall confine ourselves to briefly 
mention some papers relevant for the present work. The electric screening masses 
squared of SQED and QED have been calculated by Blaizot et al. to order and 
e^, respectively. The plasmon mass in (yf^<l>^-theory has been calculated by Parwani to 
order g'^ in resummed perturbation theory 0. The pressure has also been computed 
to order g"^ by Frenkel, Saa and Taylor 0, and extended to order g^ by Parwani 
and Singh ^j. Parwani and Coriano have carried out corresponding calculations in 
QED to order f^, and the calculational frontier has again be pushed to order 
by Parwani [|T0|. In QCD, Arnold and Zhai have obtained the free energy to order 
g^ |jTl|, while Zhai and Kastening have performed a (^^-calculation . 



Another major achievement is the application of effective field theory methods 
to the calculations of static quantities. The general idea is to take advantage of two 
or more well separated mass scales in the problem at hand, and construct a sequence 
of effective field theories which are valid at successively longer distances. The short- 
distance effects are then encoded in the parameters of the effective Lagrangian [13-15]. 
In the case of (7^$^-theory we have the scale T which is associated with the nonzero 
Matsubara frequencies (their masses are 27rnT, where n is a positive integer) and 
the scale gT which is associated with the zero-frequency mode (this mode acquires 
a thermal mass of order gT. Alternatively, this is the scale at which static scalar 
fields are screened). The strategy is to integrating out the scale T, or equivalently 
to integrating out the heavy modes, to obtain an effective Lagrangian of the light 
mode. The effective theory is three dimensional and the process is called dimensional 
reduction. In Non-Abelian gauge theories there are three momentum scales, namely 
T which again is associated with the nonstatic modes, gT which is the scale of colour 
electric screening and the scale g^T which is the scale of colour magnetic screening. 
In this case it proves useful to construct a second effective field theory by integrating 
out the scale gT 



The methods of effective field theory were first applied to high temperature field 
theory by Ginsparg and Landsman in Refs. [17,18]. In these papers, effective three- 
dimensional Lagrangians were constructed by explicitly integrating out heavy modes 
at the one-loop level. Later, Kajantie et al. have extended this method beyond the 
one-loop level by matching Greens functions in the full and in the effective theory, 
using the effective potential \T^\. This approach has mainly be used in connection 



with the study of phase transitions in the standard model and other theories [19,20]. 
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There is a nice alternative to explicitly distinguishing between light and heavy 
modes. Instead, one writes down the most general effective three-dimensional La- 
grangian which is consistent with the symmetries at high temperature. The pa- 
rameters in the effective theory are then determined by the requirement that static 
correlators in the full theory are reproduced to some desired accuracy by the effective 
theory, at long distances R ^ 1/T. This approach has been developed by Braaten 
and Nieto in Ref. [21], where they applied it to (7^$^-theory. Combining this with 
renormalization group methods, they computed the screening mass squared to order 
g^lng and the pressure to order g^lng. Later, Braaten and Nieto computed the 



free energy in QCD, through order g^ [^, and confirmed the resummation results of 



Zhai and Kastening ||T^ . The method has also been applied by the present author to 
calculate the screening mass in SQED and QED to order and e^, respectively, as 
well as the free energy in QED to order [^. These calculations have reproduced 
the results first obtained by resummation [5,9,10]. 

In this letter we apply these ideas to scalar theory and calculate the screening 
mass squared to order g^. Using the renormalization group techniques of Ref. [^, we 
also sum up leading logarithms of T/ [gT) from higher order of perturbation theory. In 
section two, we briefly discuss the effective three-dimensional theory and we determine 
the coefficients in the effective Lagrangian. In section three, we calculate the screening 
mass to order g^. In section four we summarize our results and in the appendices 
A and B, the notation and conventions are given. We also list the sum-integrals of 
the full theory and the integrals of the effective theory, which are necessary for the 
present calculations. 



2 The Effective Lagrangian 

In this section we discuss the construction of the effective three-dimensional La- 
grangian and we determine the parameters to the accuracy necessary for calculating 
the screening mass squared to order g^. Since the results in this section appear else- 
where in the literature [18,21], the discussion will be rather brief. In particular, we 
refer to the paper by Braaten and Nieto for a thorough discussion of this ap- 
proach and calculational details. 

If we denote the scalar field in the full theory by $, the Euclidean Lagrangian for 
massless ^f^^^-theory is 

CE = lid,^r+f^^'. (1) 
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In the effective theory, we similarly denote the scalar field by 0. The field in the 
three-dimensional Lagrangian can be approximately identified with the light mode in 
the underlying theory. At leading order in g"^, the relation between them is 

' (2) 



First, we must identify the symmetries of the effective Lagrangian. We have a Z2- 
symmetry — —0, which follows from the corresponding symmetry in the full theory. 
There is also a three-dimensional rotational symmetry. Hence, we may write 

£eff = \{d^ct>f + + + + /^2(A)0^VV^ + SC. (3) 

Here, we have indicated that the parameters in the effective Lagrangian generally 
depend on the ultraviolet cutoff A in the effective theory[|. This is necessary in order 
to cancel the A-dependence that arises in the calculations using the three-dimensional 
effective Lagrangian. According to Ref. [^], A can be viewed as an arbitrary factor- 
ization scale, which separates the scales T and gT. Moreover, 5C includes all other 
operators which are consistent with the symmetries. The parameters m?{K) and A(A) 
are the only operators which contribute to the screening mass squared through order 
g^. Thus for the present calculation, these are the only quantities which must be 
determined. 

The parameters in the effective theory are determined by a matching requirement. 
We demand that static correlators in the full theory be the same as those in the 
effective theory to some desired accuracy at long distances R ^ 1/T. We shall 
carry out the matching using strict perturbation theory [^. Strict perturbation 



corresponds to the following partition of the Lagrangian into a free piece and an 
interacting part: 

^0 = \{d,^)\ (4) 

Ant = Ya^'- (5) 
The Lagrangian of the effective theory is split likewise: 

(AfF)o = \{d4f. (6) 
(^eff)i,t = ^m2(A)02 + ^0^ (7) 

It is clear from Eqs. (^) and (^) that strict perturbation theory is ordinary perturba- 
tion theory, and it is therefore afflicted with severe infrared divergences, since we have 

^The field also depends on the scale A, but it is suppressed for notational ease. 
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not rearranged our Lagrangian in order to screen them. Nevertheless, this does not 
prevent us from determining the parameters using this partition of the Lagrangians. 
The point is that the parameters in C^q encode the physics on the scale T and are 
insensitive to the scale gT. So as long as we make the same incorrect assumptions in 
the effective theory, the divergences cancel in the matching procedure and we can use 
strict perturbation theory to compute the coefficients in CqQ. This implies that we are 
free to use any suitable infrared cutoff. In the present work we use dimensional reg- 
ularization. Finally, since strict perturbation theory is ordinary perturbation theory, 
it is clear that the parameters in the effective Lagrangian can be written as powers 
series in g"^. 

When we match static Greens functions in the two theories, it is necessary to take 
into account the different normalizations of the fields. The simple relation Eq. 
breaks down beyond leading order in g"^, and this is due to the wave function renor- 
malization of $ in the full theory. In gf^^^-theory, this is a two-loop effect and so is 
proportional to g^. It is therefore relevant for calculations first at order g^, and so 
we can use Eq. (^) as it stands. 

The mass parameter m^(A) must be determined at next-to- leading order in g'^ 
Denoting the self-energy function of $ by S(fco, k), we may write the static two-point 
function in the full theory as 

r(2)(A;o = 0,k) = A;2 + S(A;o = 0,k). (8) 

Here, A; = |k|. In Figs. |I] and ^ we have displayed the one and two- loop diagrams 
contributing to the self-energy function. The setting sun graph is dependent on the 
external momentum (/cq, k), but it is consistent to evaluate this diagram at vanishing 
external momentum: 

E(l, = 0, k = 0) = i - 'jg^ ^ - (9) 

Here, Zg2 is the renormalization constant for the coupling g"^. The second term 
in Eq. (|^) has a linear infrared divergence for qq = 0, while the third term has a 
logarithmic infrared divergence for po = qo = 0. Both singularities are of course 
canceled by corresponding IR-divergences in the effective theory. 

Let us now turn to the self-energy function of the field (p, which we denote by 
Sgff(k, A). It has a perturbative expansion which is given by the diagrams in Figs. |^ 
and ^ as well as the diagrams with a mass insertion. These are depicted in Fig |^. 
The self-energy function of the field is also evaluated at k = 0. This implies 
that Sgff (0, A) is identically zero, since there is no mass scale in the corresponding 
integrals. The two-point function in the effective theory can then be written as 

r|^|(k, A) = k^ + m\A) + 5m\ (10) 
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Here, we have included a mass renormalization counterterm Srn^. The mass param- 
eter m^(A) is then found by demanding 

r(2)(A:o = 0,k) = r^'f)(k,A). (11) 
This imphes that the mass parameter is given by 

m\A) = » 1,-14'-^ ' , - (12) 



The necessary sum-integrals have been listed in Appendix A. After renormalization 
of the coupling constant in Eq. (^), which is carried out by the substitution 

we are still left with a pole in e. This pole is canceled by 5m? which reads 

The mass parameter then becomes 



The mass parameter m^(A) was first obtained by Braaten and Nieto in Ref. PT|. We 
have used the renormalization group equation for the coupling g'^ 

^ = 16^' 

to change the scale from A to and so the coupling constant g is then evaluated at 
the scale \i. 

It is clear from Eq. (|I5|) that m^(A) depends explicitly on the arbitrary scale A. 
The mass parameter satisfies an evolution equation, which can be derived from the 
requirement that the physical screening mass be independent of this arbitrary scale 
A |21|. Alternatively, one can differentiate Eq. (^) with respect to A. In terms of 
A (A) the evolution equation reads 

'^-iK = 6(4^) • 

The physical interpretation of the mass parameter is that it is contribution from 
the scale T to the screening mass. 
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Let us now move on to the quartic self-interaction. We need to know the coeffi- 
cient A (A) in front at next-to-leading order in g"^. Using similar arguments as we did 
when we computed the mass parameter, one finds 



\{K) = Zg.g^T--g'T 



1 

P P4- 



The first term on the right hand side results from matching at tree level, while the 
second term arises from matching at the one-loop level. The corresponding diagram 
is displayed in Fig. H. After renormalization of g"^, we find 



A(A)=/T- 



?,g^T 



(477) 



In 



A 



+ 7e 



(19) 



This coefficient was ffist found by Landsman in Ref. |T8[ by explicitly integrating out 
the nonzero Matsubara modes. Using the renormalization group equation for g"^ one 
can demonstrate that the quartic coupling in the effective theory is independent of 
the scale A, and so we can trade it for an arbitrary scale fi. 



3 The Screening Mass 

In this section we calculate the screening mass squared to order g^ using the effective 
Lagrangian. When we calculated the parameters m^(A) and A(A), we treated the in- 
frared divergences incorrectly by decomposing our Lagrangians according to Eqs. (^ 
and (H), and Eqs. (||) and (0). In order to treat these divergences correctly, we must 
now incorporate the infrared cutoff provided by the mass parameter into the free part 
of £eff. Thus, we split the effective Lagrangian according to 

(Aff)o = l^d^f + lm'm^ (20) 

(^eff)int = ^<P'- (21) 

The screening mass is defined as the pole of the propagator at spacelike momen- 
tum [grj: 

ml = m\A) + Scff(k, A), = -m^. (22) 

The Feynman diagrams contributing to the self-energy function Sgg-(k, A) at next- 
to-next-to-leading order are depicted in Figs. H — and |^: 

Seff(k,A) = - / — --6m^- 



2 jp p2 _|_ ^2 2 Jp {p^ + m^Y 4 Jpq {p^ + m^Y^q^ + 



6 
8 

A3 
8 

A3 
8 

A3 
12 



1 



pq {p^ + m?) {q^ + m?) [(p + q + k)^ + m?] 
1 1 1 



+ 



pqr 



pqr 



1 



+ 



1 



+ 



pgr (p^ + mP'Y^q^ + m^)[(p + q + k)^ + m2](r^ + m^^ 

1 



+ 



pqr 



6m^. (23) 



The mass counterterm cancels the pole in e from the second two-loop diagram 
in Fig. 1^. Similarly, the divergence of the one-loop diagram with a mass counterterm 
insertion cancels against a pole in e arising from the third three-loop graph of Fig. |^. 
Moreover, this one-loop diagram also gives finite contributions to the screening mass, 
when 1/e "hits" the terms proportional to e. Finally, the last two-loop diagram 
and the last three-loop diagram depend on the external momentum k. The effective 
expansion parameter of the quartic interaction is A/m, and in order to compute 
the screening mass squared consistently in powers of A/m, we must evaluate these 
integrals at = im. In order to see that this in fact is necessary, one can perform 
an expansion in the external momentum k, and verify that all terms are equally 
important for soft k ~ gT 

The integrals appearing in the expression for the screening mass are tabulated in 
Appendix B. Expanding the mass parameter in powers of g"^, we obtain the screening 
mass squared through order g^: 



2 rji2 J 1 

m: = — 1 < 1 
24 1 



.C'(-i) 



C(-i 
C'(-i) 



+ 



47r (47r)2 

.3 .9 



31n-^ + 41n- 



AtiT 



AttT 



1 + 81n2-7ij + 2 + 
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(47r) 



:V6 



u , 2m 
- In — 2 In 

2 AnT AnT 



3-71n2 + -7i5 



C(-i) 



(24) 



Firstly, one notes that the A-dependence explicitly cancels in Eq. ( ^4]) to next-to- 
next-to-leading order in g. Using the running of the coupling g"^, it also easy to verify 
that there is no dependence on the renormalization scale either. The result to order 
g^ was first obtained by Dolan and Jackiw |l[], while Braaten and Nieto computed 
the screening mass squared to order g'^. Moreover, they used the evolution equations, 
which the parameters in the effective theory satisfy, to sum up leading logarithms of 
higher order in the perturbation expansion. The result includes in particular a term 
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proportional to g^lng. This term is of course also present in Eq. p^ ) above. The 
complete result to order is new. 



We would also like to point out that the effective field theory approach explains 
the very appearance of these logarithms. The mass parameter, and possibly other 
operators as well, depend explictly on the factorization scale A and such terms occur 
as logarithms of A/T. In the effective theory, logarithms of m/A arise in perturba- 
tive calculations. In order to cancel the A-dependence in physical qauntities, these 
logarithms must match, leaving logarithms of T/m or T/{gT). Hence, these loga- 
rithms can be attributed to the renormalization of the parameters in the effective 
Lagrangian , and in the present case it is the renormalization of the mass param- 



eter. 

The evolution eqution for the mass parameter, Eq. (0), can be used to sum up 
the leading logarithms of T/ (gT). This can be done by choosing the renormalization 
scale fj, = 2ttT and the factorization scale A = in the expression for the mass 
parameter The screening mass squared again follows from Eqs. (|2^) and (p3D , 

and in terms of and A it reads 

ml = m2|l-l(^)-(-^)2[i_lln2]-(-^f[- + -ln2]|. (25) 

This is the full result to order g^ and also includes leading logarithms in the form 
^2n+3|j^n^^ where n is a natural number. These terms are obtained by expanding the 
second term in Eq. (E5|) in powers of g^. 



4 Summary 

In the present work we have obtained the screening mass squared to order g^ using 
effective field theory methods. The short-distance coefficients contains the physics on 
the scale T, while the effective theory takes care of contributions to physical quantities 
from the scale gT. Thus, effective field theory methods unravel the contributions to 
physical quantities from the different momentum scales, and streamlines calculations, 
since we treat one scale in the problem at a time. This is the advantage of the 
effective field approach over the more conventional resummation procedure; the latter 
complicates calculations unnecessarily because the sum-integrals involve both scales 
T and m. 

We were also able to sum up leading logarithms of T/{gT) arising from higher 
orders of perturbation theory by using the evolution equation satisfied by the mass 
parameter. 
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Our result is correct up to corrections of order . In order to obtain the screening 
mass squared to this order, we must determine the mass parameter to three-loop 
order and also compute the four-loop diagrams in the effective theory. There are 
no new coefficients in the effective Lagrangian which must be determined, since the 
corresponding operators first contribute to at order or higher. The calculation 
should be a manageable task using the machinery developed to handle difficult multi- 
loop sum- integrals lO, and complicated integrals in the effective theory. 



A Sum- integrals in the Full Theory 



In this appendix we give expressions for the sum-integrals used in the present work. 
We use the imaginary time formalism, where the four-momentum is P = (po? p) 
with = Po + P^- "^^^ Euclidean energy takes on discrete values, po = '^^unT for 
bosons. Dimensional regularization is used to regularize both infrared and ultraviolet 
divergences by working in c? = 4— 2e dimensions, and we apply the MS renormalization 
scheme. We use the following shorthand notation for the sum-integrals that appear 
below: 



47r 



P0=27rnT 



(27r; 



Then /i coincides with the renormalization scale in the MS renormalization scheme. 

Arnold and Zhai have developed the machinery necessary to evaluate complicated 
multi-loop sum-integrals [Tl|. They have calculated and listed the specific sum- 
integrals needed in the present work, and details may be found in Ref. lTl|]. We 
list them here for the convenience of the reader 



1 

PP2 



12 



1 + (2 In 



p {P 



2\2 



rpQ p'^Q'^iP + qy 



(4vr) 
0. 



AtiT 



2 + 2 



- + 21n 

e AttT 



C(-i 

+ 2^E + 0{e 



(A.2) 
(A.3) 
(A.4) 



B Integrals in the Effective Theory 

In the effective three-dimensional theory we use dimensional regularization in 3 — 
2e dimensions to regularize infrared and ultraviolet divergences. In analogy with 
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Appendix A, we define 



P 



3~2e 



fip)- 



(B.i; 



Again /i coincides with the renormahzation scale in the modified minimal subtraction 
renormalization scheme. All the integrals necessary in the present work, except for 
those in Eqs. ([B.4[ ), (|B.6| ) and ( [B.7| ) below, have been calculated by Braaten and 
Nieto in Ref. I have computed these remaining integrals using their methods. 



p2 _|_ ^2 



1 



m 
Air 
1 



l + (21n-^ + 2)e + 0(e2) 
l+(21n-^)e + 0(e^ 



1 + (21n^ + 2)e 

V ).m ' 



2m 



L 



pq {p^ + m^) (g2 + m^) [(p + q + k)^ + m 



k=im 

^ 'i + 6 + 41n-^-81n2 + 0(e) 
e 2m 



in 



1 



pg {p^ + m^Y^q^ + m^) [(p + q + k)^ + m^ 

1 



1 



k=im 

ln2 + 0(e 



pqr (p^ + m'^Y{q'^ + m'^){r'^ + m'^)[{p + q + r)^ + 



1 



m(87r)' 



i + 61n — - 41n2 + 2 + 0(e) 
e 2m 



(B.2) 
(B.3) 
(B.4) 



(B.5) 



(B.6) 



(B.7) 



References 

[1] L. Dolan and R. Jackiw, Phys. Rev. D 9, 3320, 1974. 

[2] E. Braaten and R. Pisarski, Nucl. Phys. B 337, 569, 1990. 

[3] M. Le Bellac, Thermal Field Theory, Cambridge University Press, Cambridge 
1996. 

[4] F. Flechsig and A. K. Rebhan, Nucl. Phys. B 464, 279, 1996. 

[5] J.-P. Blaizot, E. lancu and R. Parwani, Phys. Rev. D 52, 2543, 1995. 



11 



[6] R. Parwani, Phys. Rev. D 45, 4695, 1992. 

[7] J. Frenkel, A. V. Saa and J. C. Taylor, Phys. Rev. D 46, 3670, 1992. 

R. Parwani and H. Singh, Phys. Rev. D 51, 4518, 1995. 

[9] R. Parwani and C. Coriano, Phys. Rev. Lett. 73, 2398, 1994. 

[10] R. Parwani, Phys. Lett. B 334, 420, 1994, Erratum-ibid. B 342, 454, 1995. 

[11] P. Arnold, C.-X. Zhai, Phys. Rev. D 50, 7603, 1994; Phys. Rev. D 51, 1906, 
1995. 

[12] C.-X. Zhai and B. Kastening, Phys. Rev. D 52, 7232, 1995. 

[13] D. Kaplan, |nucl-th / 9506035i Lectures given at 7th Summer School in Nuclear 
Physics Symmetries. 



[14] A. Manohar, |hep-ph /9508245| . Lectures given at 1995 Lake Louise Winter Insti- 



tute on Quarks and Colliders. 

[15] G. P. Lepage, in From actions to answers, proceedings of the Advanced Study 
Institute in Elementary Particle Physics, Boulder, Colorado, 1989, edited by T. 
De-Grand and D. Toussaint, World Scientific, 1989. 

[16] E. Braaten, Phys. Rev. Lett. 74, 2164, 1995. 

[17] P. Ginsparg, Nucl. Phys. B 170, 388, 1980. 

[18] N. P. Landsman, Nucl. Phys. B 322, 498, 1989. 

[19] K. Kajantie, M. Laine, K. Rummukainen and M. Shaposhnikov, Nucl. Phys. B 
458, 90, 1996. 

[20] D. Bodeker, P. John, M. Laine, and M. G. Schmidt, Nucl. Phys. B 497, 387, 
1997. 

[21] E. Braaten and A. Nieto, Phys. Rev. D 51, 6990, 1995. 
[22] E. Braaten and A. Nieto, Phys. Rev. D 53, 3421, 1996. 
[23] J. O. Andersen, Phys. Rev. D 53, 7286, 1996; Z. Phys. C 75, 147, 1997. 



12 



Q 



Figure 1: One-loop correction to the two-point function. 



Figure 2: Two-loop corrections to the two-point function. 




Figure 3: Diagrams with a mass insertion contributing to the two-point function in 
the effective theory. 
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Figure 4: One-loop correction to the four-point function. 




Figure 5: Three-loop corrections to the two-point function in the effective theory. 
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